Preprint typeset in JHBP style - HYPER VERSION 



Reconciling thermal leptogenesis with the gravitino 
problem in SUSY models with mixed axion/axino dark 
matter 



Howard Baer a , Sabine Kraml 6 , Andre Lessa a and Sezen Sekmen c 

a Dept. of Physics and Astronomy, University of Oklahoma, Norman, OK 73019, USA 
b Laboratoire de Physique Subatomique et de Cosmologie, UJF Grenoble 1, 
CNRS/IN2P3, INPG, 53 Avenue des Martyrs, F-38026 Grenoble, France 
c Dept. of Physics, Florida State University, Tallahassee, FL 32306, USA 



E-mail: baerOnhn . ou . edu 



sezen . sekmenOcern . ch 



sabine . kramlQlpsc . in2p3 . f r , lessa. a.p@gmail . com 



Abstract: Successful implementation of thermal leptogenesis requires re- heat tempera- 
tures Tr >2x 10 9 GeV, in apparent conflict with SUSY models with TeV-scale gravitinos, 
which require much lower Tr in order to avoid Big Bang Nucleosynthesis (BBN) constraints. 
We show that mixed axion/axino dark matter can reconcile thermal leptogenesis with the 
gravitino problem in models with > 30 TeV, a rather high Peccei-Quinn breaking 
scale and an initial mis-alignment angle 9i < 1. We calculate axion and axino dark matter 
production from four sources, and impose BBN constraints on long-lived gravitinos and 
neutralinos. Moreover, we discuss several SUSY models which naturally have gravitino 
masses of the order of tens of TeV. We find a reconciliation difficult in Yukawa-unified 
SUSY and in AMSB with a wino-like lightest neutralino. However, Tr ~ 10 10 — 10 12 GeV 
can easily be achieved in effective SUSY and in models based on mixed moduli-anomaly 
mediation. Consequences of this scenario include: 1. an LHC SUSY discovery should be 
consistent with SUSY models with a large gravitino mass, 2. an apparent neutralino relic 
abundance ^^^- 2 1, 3. no WIMP direct or indirect detection signals should be found, 
and 4. the axion mass should be less than ~ 10~ 6 eV, somewhat below the conventional 
range which is explored by microwave cavity axion detection experiments. 



Keywords: 
Matter. 



Supersymmetry Phenomenology, Supersymmetric Standard Model, Dark 



1. Introduction 



Recent measurements of neutrino oscillations M are elegantly interpreted in terms of see- 
saw neutrino masses ||], wherein heavy right-handed neutrino (RHN) states Ni (i = 1—3 for 
three generations) are introduced, and the light neutrino masses are given approximately 
by m Ui ~ m^./MjVj, where mm ~ f Vi v with f„. the neutrino Yukawa coupling and v the 
Higgs field vacuum expectation value. A value of ~ 10 15 GeV yields m Vr ~ 0.03 eV 
in the GUT-inspired case where f Ur = ft at Mqut- 

One of the appealing consequences of such heavy RHN states is that the baryon number 
of the universe can be explained in terms of thermal leptogenesis ||. In thermal leptoge- 
nesis, the right-hand neutrinos are present in thermal equilibrium at high temperatures in 
the early universe, and decay asymmetrically to leptons versus anti-leptons. The lepton 
asymmetry is converted to a baryon asymmetry via B- and L-violating but B — L con- 
serving sphaleron interactions. In a scheme with hierarchical right-hand neutrino masses, 
a re-heat temperature 

T R > 2 x 10 9 GeV (thermal leptogenesis) (1.1) 

is required to reproduce the observed baryon asymmetry of the universe @. 

The existence of a heavy mass scale like Mn or Mqut brings about the infamous 
hierarchy problem of the Standard Model (SM). The gauge hierarchy problem is elegantly 
resolved by introducing supersymmetry (SUSY) into the theory. SUSY is a novel spacetime 
symmetry which relates bosons and fermions, see e.g. ||. In realistic models, SUSY is 
broken "softly" at the weak scale — implying that all SM particles must have superpartners 
with masses in the range up to 0(1) TeV that ought to be accessible to colliders such 
as the CERN LHC. SUSY not only stabilizes the hierachy between weak scale and other 
high scales such as Mqut; indeed with weak-scale SUSY the gauge couplings, when evolved 
upwards from Q = Mz under renormalization group evolution (RGE) , are found to unify at 
Q ~ 2x 10 16 GeV, which is nicely consistent with the idea of a grand unified theory. Sensible 
implementations of supersymmetry invoke SUSY as a local symmetry (supergravity or 
SUGRA), which requires in addition the existence of a graviton/gravitino supermultiplet. 
In SUGRA models, supergravity is broken via the superHiggs mechanism, leading to a 
massive gravitino G. The soft SUSY breaking mass terms are expected to be closely 
related to the gravitino mass mg, and so mg is also expected to be of the order of the 
weak scale. 

One of the impediments to successful SUGRA model building is known as the gravitino 
problem Q. Gravitinos can be produced thermally in the early universe, even though their 
Planck-suppressed couplings preclude them from participating in thermal equilibrium. The 
gravitino decay rate is also suppressed by the Planck scale, leading to very long gravitino 
lifetimes of order 1 — 10 5 sec. There are actually two parts to the gravitino problem. 
Part 1 is that for ~ 1 TeV, the late time gravitino decays can inject hadronic or 
electromagnetic energy into the cosmic plasma at a time scale during or after Big Bang 
Nucleosynthesis (BBN), leading to destruction of the successful agreement between theory 
and observation for the light element abundances. For re-heat temperatures Tr < 10 5 GeV, 
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thermal gravitino production is suppressed enough to evade BBN limits @. However, these 
low of temperatures are in conflict with those needed for thermal leptogenesis. 

If 5 TeV < < 30 TeV, then the gravitino lifetime drops below 1 sec, and the BBN 
constraints are much more mild, allowing T R < 10 9 GeV. This range of Tr is consistent 
with non-thermal leptogenesis ||, wherein other sources of N{, such as inflaton decay, 
contribute to Ni production. For heavier yet gravitinos with mg > 30 TeV, the value of Tr 
can reach as high as 7 x 10 9 GeV. In this case — part 2 of the gravitino problem — the upper 
bound on Tr comes from overproduction of neutralino dark matter due to their combined 
thermal production and production via gravitino decay. Models such as AMSB, with 
multi-TeV gravitino masses and very low thermal neutralino abundances, thus naturally 
reconcile thermal leptogenesis with the gravitino problem, but only for the narrow range 
2 x 10 9 GeV<TR< 7x 10 9 GeV §. 

One way out is to invoke a gravitino as LSP, with a stau or neutralino as the next-to- 
LSP (NLSP) which decays via a small i?-parity violating interaction [jl(| (in the R-parity 
conserving case, it is very difficult to reconcile gravitino DM with thermal leptogenesis [10 



[Tl[]). The gravitino, which may also decay via i?-parity violating interactions, has a lifetime 
of order the age of the universe; it can still function as dark matter, but its occasional 
decays in the galactic halo could be the source of Pamela, ATIC and Fermi cosmic ray 
anomalies Jl^]. 

Another way out involves mixed axion/axino DM, and this is the topic of this paper. 
Indeed, the strong CP problem remains as one of the central puzzles of QCD which evades 
explanation within the context of the Standard Model. The crux of the problem is that an 
additional CP violating term in the QCD Lagrangian of the form 1 9gl / '32-7T 2 G A ought 
to be present as a result of the t'Hooft resolution of the U{\)a problem via instantons and 



the 9 vacuum of QCD [13]. The experimental limits on the neutron electric dipole moment 
however constrain \9\ < 10~ 10 113] . Why this term should be so small is the essence of the 
strong CP problem. 



An extremely compelling solution proposed by Peccei and Quinn [15] is to hypothesize 
an additional global ?7(1)pq symmetry, which is broken at some high mass scale f a . A 
consequence of the broken PQ symmetry is the existence of a pseudo-Goldstone boson 
field: the axion a(x) [l(J. The Lagrangian then also contains the terms 



C ^2 d ^ a+ ^m G ^- <L2> 



where ./V is the model-dependent color anomaly factor, which is 1 for KSVZ [17] or 6 



for DFSZ [18] models. Since a(x) is dynamical, the entire CP- violating term settles to 
its minimum at zero, thus resolving the strong CP problem. A consequence of this very 
elegant mechanism is that a physical axion field should exist, with axion excitations of 
mass [19] 

ma ~ 6 eV • 



1 Here G„„ is the gluon field strength tensor and Q a ^ v its dual. 
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The axion field couples to gluon-gluon (obvious from Eq. ( |1.2[) ) and also to photon-photon 
and fermion-fermion. All the couplings are suppressed by the PQ scale f a . The value of f a 
is constrained to lie above ~ 10 9 GeV by stellar cooling arguments pCfl , leading to a nearly 
invisible axion particle which may be searched for via microwave cavity experiments plfl . 
In addition, axions can be produced via various mechanisms in the early universe. Since 
their lifetime (they decay via a — > 77) turns out to be longer than the age of the universe, 
they can be a good candidate for the DM of the universe 22]. 



In the context of supersymmetry, the axion field is but one element of an axion su- 
permultiplet which also contains an i?-parity even spin-0 saxion field s(x) and an .R-parity 
odd spin-^ axino field a [23]. The axino field a may play a huge role in cosmology [24]: its 



mass may lie anywhere in the range of keV to TeV 25], and it may function as the LSP. As 



the LSP, in i?-parity conserving models, it may constitute at least a portion of the DM of 



the universe [26, 27]. The saxion field may also play a role in cosmology, e.g. via dilution 



of relics by additional entropy production [28, E9], although we will not consider this here. 



In [30], Asaka and Yanagida proposed to reconcile thermal leptogenesis with the grav- 
itino problem by requiring an axino LSP with a gravitino NLSP. In this paper, we present 
a different reconciliation of thermal leptogenesis with the gravitino problem, nota bene for 
very heavy gravitinos. First, we require > 30 TeV so as to avoid part 1 of the gravitino 
problem. Next, we invoke the presence of mixed axion/axino dark matter into our scenario, 
which arises as a result of the PQ solution to the strong CP problem in the supersymmet- 



ric context [31|. We assume here that the lightest neutralino, Z\, is the NLSP, so that 
each neutralino ultimately decays to an axino, and the neutralino relic mass abundance is 
then suppressed by a factor of mj/m^, thus avoiding part 2 of the gravitino problem, the 
overproduction of neutralino dark matter. 

In Sec. H of this paper, we evaluate the relic abundance of mixed axion/axino DM due 
to four sources. The first is ordinary production of axion cold DM via the vacuum mis- 
alignment mechanism. The second is thermal production (TP) of axinos: here, we restrict 
Tr to values below the axino decoupling temperature (T^ cp ), so axinos are never in thermal 
equilibrium. Nonetheless, axinos can still be produced thermally via bremsstrahlung and 
decays of particles which are in thermal equilibrium: the final result depends linearly on 
the re-heat temperature Tr after inflation. The third is non-thermal production (NTP) 
of axinos via thermal neutralino production and decay. The fourth, also NTP of axinos 
comes from thermal production of gravitinos, followed by their cascade decays to the axino 
LSP state; this mechanism also depends linearly on Tr. 

For low values of PQ breaking scale f a /N ~ 10 9 — 10 11 , the axino coupling to matter is 
large enough that thermal production of axinos tends to dominate the mixed axion/axino 
abundance. An upper bound on Tr can be extracted by requiring the axino DM abundance 



lie below the WMAP-measured value [32|: 

n DM h 2 = 0.1123 ± 0.0035 at 68% CL. (1.4) 

As fa/N increases, the portion of axino DM decreases, while the axion CDM increases. 
Simplistic estimates of the relic axion abundance assume an initial mis-alignment angle 
9i ~ 1, leading to an upper bound on PQ breaking scale f a /N < 5 x 10 11 GeV. By 
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adopting a smaller value of 0j, much larger values of f a /N in the 10 13 10 14 GeV range 
become allowed. This in turn suppresses the axino composition of DM, unless very high 
values of Tr > 10 10 GeV are allowed. This is the crux of our reconcilation of thermal 
leptogenesis with the gravitino problem. 

However, within this solution, another BBN bound emerges, since the lifetime of the 
lightest neutralino scales as (f a /N)~ 2 . There exist additional strict limits on late decaying 
particles with electromagnetic and hadronic energy injection into the thermal plasma during 
BBN. Therefore, these limits provide an upper bound on f a /N for SUSY models including 
the PQ mechanism. In Sec. ||, we evaluate the neutralino lifetime and hadronic branching 
fraction, so that BBN constraints can be applied to Z\ decay. 

In Sec. ||, we present five scenarios which are consistent with gravitinos of mass > 30 
TeV. The first two cases come from gravity mediated SUSY breaking: the Yukawa-unified 
(YU) models and effective SUSY (ESUSY) models: both of these require GUT scale scalar 
masses in the 10 — 30 TeV regime, and so should be consistent with gravitinos in this 
range. As we will see, the rather large abundance of bino-like neutralinos from YU models 
typically makes them inconsistent with Tr values in excess of 2 x 10 9 GeV. ESUSY models 
can more easily accommodate a low production rate for neutralinos in the early universe, 
and do allow for T R > 10 10 GeV. The third case, that of anomaly-mediated SUSY breaking 
(AMSB), requires gravitinos in the 30 — 100 TeV range since soft SUSY breaking terms 
are loop suppressed. Most versions of AMSB include a nearly pure wino-like neutralino. 
Since only the bino component of Z\ couples to the axino, the neutral wino-like Z\ decay 
is suppressed, and just barely allows Tr > 2 x 10 9 GeV, before BBN constraints kick in. 
The fourth and fifth benchmark points come from "mirage unification" (MU), or mixed 
moduli-anomaly mediation, which also allows ~ 30 — 100 TeV. These models can easily 
include a bino-like Z\ with a relatively low relic abundance, avoiding the worst of BBN 
constraints. In these models, values of Tr ~ 10 9 — 10 12 can easily be accommodated, 
thus reconciling thermal leptogenesis with the gravitino problem. In Sec. [5L we present a 
summary and conclusions. 

2. Mixed axion/axino relic density 

In this section, we list the four production mechanisms for mixed axion/axino dark matter 
which are considered here. It is possible that other more exotic mechanisms could also 
contribute, such as axino production from moduli or inflaton decay. We will not consider 
these additional mechanisms here. 

2.1 Axions via vacuum misalignment 

Here, we consider the scenario where the PQ symmetry breaks before the end of inflation, so 
that a nearly uniform value of the axion field 9i = a(x)/(f a /N) is expected throughout the 
universe. The axion field equation of motion implies that the axion field stays relatively 
constant until temperatures approach the QCD scale Tqcd ~ 1 GeV. At this point, a 
temperature-dependent axion mass term turns on, and a potential is induced for the axion 
field. The axion field rolls towards its minimum and oscillates, filling the universe with 
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low energy (cold) axions. The expected axion relic density via this vacuum mis-alignment 
mechanism is given by [33] 



where < 9{ < ir and is the so-called anharmonicity factor. Visinelli and Gondolo [33 

parametrize the latter as = In ^ i-g±/ n 2 ) j • The uncertainty in £l a h 2 from vacuum 

mis-alignment is estimated as plus-or-minus a factor of three. 

2.2 Thermal production of axinos 

If the reheat temperature Tr exceeds the axino decoupling temperature 



^"^(ira) 2 ^ 3 , (-) 



axinos will be in thermal equilibrium, with an abundance given by &£ E h 2 ~ 2l?V" ^° 



avoid overproducing axino dark matter, the RTW bound [24| then implies that m„ < 0.2 
keV. We will here consider only reheat temperatures below T^ cp . In this case, the axinos 
are never in thermal equilibrium in the early universe. However, they can still be produced 



thermally via radiation off of particles that are in thermal equilibrium |26], 34]. Here, we 



adopt a recent calculation of the thermally produced axino abundance from Strumia [35|: 

,2.3, 

with F(gz) ~ 20g| In J^-, and g% is the strong coupling constant evaluated at Q = Tr. 

2.3 Axinos via neutralino decay 

In supersymmetric scenarios with the neutralino as a quasi-stable NLSP, the Z\s will be 
present in thermal equilibrium in the early universe, and will freeze out when the expansion 
rate exceeds their interaction rate, at a temperature roughly Tf ~ /20. The present 
day abundance can be evaluated by integrating the Boltzmann equation. Several computer 
codes are available for this computation. Here we use the code IsaReD |36|], a part of the 
Isajet/Isatools package [37, 38 1. 

In our case, each neutralino will undergo decay to the stable axino LSP, via decays 
such as Z\ — > aj. Thus, these non-thermally produced axinos will inherit the thermally 
produced neutralino number density, and we will simply have [^(| 

n z h 2 = H^sF/h 2 . (2.4) 

m z 1 

2.4 Axinos from gravitino cascade decay 

Since here we are attempting to generate reheat temperatures Tr > 10 9 GeV, we must 
also include in our calculations the thermal production of gravitinos in the early universe. 
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We here follow Pradler and Steffen, who have estimated the thermal gravitino production 
abundance as |3^] 



^ - E«-f ™ (. + ^) * (jAg) (^) (^) . P ,, 

2 — 1 \ Cr / 

where = (0.018,0.044,0.117), = (1.266,1.312,1.271), g { are the gauge couplings 
evaluated at Q = Tr and Mi are the gaugino masses also evaluated at Q = Tr. For 
the temperatures we are interested in, this agrees with the calculation by Rychkov and 
Strumia[p9[| within a factor of about 2. 

Since each gravitino cascade decays ultimately to an axino LSP, the abundance of 
axinos from gravitino production is given by 

mg G 

For axino masses in the MeV range and gravitino masses in the 30 — 50 TeV range, the 
prefactor above is extremely small, and allows us to evade overproduction of dark matter 
via thermal gravitino production (part 2 of the gravitino problem) . 2 3 

2.5 Mixed axion/axino dark matter 

In this paper, we will evaluate the mixed axion/axino relic density from the above four 
sources: 

vt a - a h 2 = n a h 2 + nfh 2 + n?h 2 + nfh 2 . (2.7) 

Over much of parameter space, if is taken to be of order the MeV scale, then the 
contributions from Z\ and gravitino production are subdominant. In Fig. [j], we illustrate 
in the upper frame the relative importance of the four individual contributions as a function 
of f a /N, for an mSUGRA scenario with mg = 1 TeV, = 122 GeV, VL^h 2 = 9.6 (as in 
Ref. p0|| ). For the axion/axino sector we take 0i = 0.05 and m- a = 100 keV. The value of 
Tr is adjusted such that Q a a,h 2 = 0.1123. For low f a /N values, the TP axino contribution 
is dominant. But as f a /N increases, the axion component grows until at f a /N ~ 4 x 10 13 
GeV it becomes dominant, and for even higher f a /N it saturates the DM relic density. 

The value of Tr which is needed is shown in the lower frame of Fig. [l]. We see that Tr 
grows quickly with increasing f a /N. This is because the thermal axino production decreases 
as the inverse square of f a /N, so larger values of Tr are needed to keep Q a ah 2 = 0.1123. We 
see that Tr can reach ~ 10 11 GeV in the case of mainly axion CDM (similar to Ref. pOfl ). In 

2 We have checked that for gravitinos in the mass range 5 TeV < mg < 50 TeV, the temperature Ta ecay = 
y/TQmpi I (7T 2 (7»/90) 1/ ^ 4 at which gravitinos decay ranges between 0.01-0.5 MeV, well after neutralino freeze- 
out. This means that the gravitino cascade decays indeed contribute to the ultimate axino relic density, 
rather than having the neutralinos from the cascade reprocessed in the thermal bath, as would occur for 
T decay > m 2l /20. 

3 Along with gravitino cascade decays, the decay G — >• aa can also occur, and give rise to a component of 
hot dark matter axions. Since this contribution to the axion density is ~ (m a /mg)f2g/i 2 with m a ~ M e V, 
it will be a very tiny contribution to the total dark matter density, and so we are safe to neglect it here. 
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m_ = 100 keV, m_ = l TeV, m. = 122 GeV, Q r h 2 = 9.6, 9. = 0.05 
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Figure 1: Upper frame: Contribution of axions and TP and NTP axinos to the DM density as a 
function of the PQ breaking scale f a /N, for an mSUGRA point with mo = 1000 GeV, m 1 / 2 — 300 
GeV, Aq — 0, tan/3 = 10 and /i > 0, and fixing rria = 100 keV and B. L — 0.05; Tr is adjusted such 
that flaah 2 = 0.1123. Lower frame: the Tr that is needed to achieve £l a ah 2 = 0.1123 for = 0.1 
and 1 MeV, for the same mSUGRA point and 8i. 



our case here, allowing a smaller value of 0, allows higher values of f a /N to be found, which 
in turn requires much higher values of Tr, into the range needed for thermal leptogenesis. 

However, for such high / a /-/V, the Z\ becomes so long-lived that it violates the bounds 
from BBN on late decaying neutral particles (as indicated in the figure). We address this 
issue in the next section. 

3. Neutralino lifetime and hadronic branching fraction 

We have averted one problem with BBN by requiring the presence of a gravitino with 
m G ~ 3® TeV, so that it will decay largely before BBN starts. In the process, by asking 
for Tr > 2 x 10 9 GeV while avoiding overproduction of mixed axino/axion dark matter 
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(the latter requires large f a /N ~ 10 12 GeV and small 9i), we have pushed the Z\ lifetime 
uncomfortably high, so that its hadronic decays in the early universe now have the potential 
to disrupt BBN. 

Constraints from BBN on hadronic decays of long-lived neutral particles in the early 
universe have been calculated by several authors [41, 42, 43]. Here, we will adopt the 



results from the recent calculations by Jedamzik [43|. The BBN constraints arise due to 



injection of high energy hadronic particles into the thermal plasma during or after BBN. 
The constraints depend on three main factors: 



• The abundance of the long-lived neutral particles. In Ref. [43], this is given by Qxh 2 
where X is the long-lived neutral particle which undergoes hadronic decays. In our 
case, where the long-lived particle is the lightest neutralino which decays to an axino 
LSP, this is just given by the usual thermal neutralino abundance h 2 , as calculated 
by IsaReD @. 

• The lifetime tx of the long-lived neutral particle. Obviously, the longer-lived X is, 
the greater its potential to disrupt the successful BBN calculations. 

• The hadronic branching fraction Bh of the long-lived neutral particle. If is very small, 
then very little hadronic energy will be injected, and hence the constraints should be 
more mild. 

The BBN constraints are shown in Fig. 9 (for rax = 1 TeV) and Fig. 10 (for mx = 100 
GeV) of Ref. (|43fl , as contours in the tx vs. VLxh 2 plane, with numerous contours for 
differing Bh values ranging from 10~ 5 to 1. For Bh ~ 0.1, for instance, and very large 
values of Vixh 2 ~ 10 — 10 3 , the lifetime tx must be < 0.1 sec, or else the primordial 
abundance of 4 He is disrupted. If Vixh 2 drops below ~ 1, then much larger values of tx 
up to ~ 100 sec are allowed. If one desires a long-lived hadronically decaying particle in the 
early universe with tx > 100 sec, then typically much lower values of £lxh 2 ~ 10~ 6 — 10 -4 
are required. Such low neutralino relic densities are extremely hard to generate in SUSY 
models, even in the case of AMSB, or pure higgsino annihilation [44]. We have digitized the 



constraints of Ref. ^3|, implementing extrapolations for cases intermediate between values 
of parameters shown, so as to approximately apply the BBN constraints to our scenario 
with a long-lived neutralino decaying hadronically during BBN. 

3.1 Neutralino lifetime 

We have calculated the two-body decays of Z\ to axinos, and find agreement with the 
results of Ref. [26]. In the notation of Ref. ||, we find 



and 



V( y , a _ a2 Y C lYY ^ C ^ °W 3 (1 m l \ ,0 n 

" 128VT3 { fa/N) 2 y-^J (3 ' 1} 



r(Z t -4 ~aZ) = a v Cl Y V ^\l W rn% A§(1, ^ ml 



128vr 3 (fa/N) 2 Zi K ' m 2 ~ ' m 2 
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m,h 



2 



(1 



m 



a x2 



+ 



2m 2 ~ 
Zi J 



(3.2) 



with ay = (e 2 /4vr)/cos 2 (9^ the ?7(l)y coupling; C aYY = 8/3 in the DFSZ model @ and 
0, 2/3 or 8/3 in the KSVZ model |l7j depending on the heavy quark charge eQ = 0, —1/3 
or 2/3. Throughout our analysis we assume C a YY = 8/3. The above decays should be the 
only two-body decay modes allowed for the KSVZ model; for the DFSZ model, additional 
decays to higgs states may also be allowed. Since such decays dependent on the type of 
DFSZ model, we do not consider them in our analysis. 

Using the above formulae, in Fig. || we plot the lifetime t{Z\) in seconds versus 

for various choices of (f a /N)/v^ . The quantity v± denotes the bino-component of Z\ 
in the notation of Ref. ||. For models with a bino-like Z%, v± ~ 1. For models with a 
wino-like Z\, as in AMSB, the lifetime will be enhanced by a large factor, since in these 
models v$ is typically 10~ 2 — 1CP 3 . 

From Fig. ||, we see that for models with a bino-like neutralino and t(Z\) < 0.01 sec, 
either very small f a /N < 10 10 GeV are required, or if larger f a /N values are desired, then 
nig must be very (perhaps uncomfortably) large. However, if t(Z\) < 10 2 sec is needed, 
then values of f a /N as large as 10 14 GeV are allowed, depending on . 



3.2 Two- and three-body Z\ decay to hadrons 

Finally, to implement BBN constraints, we will need the hadronic branching fraction of Z\ 
decay. If the decay Z\ — > dZ is open, then Bh is just given by Bh = v ^ z y^ aZ '> x BF(Z — > 

hadrons). When the decay Z\ — > dZ is closed, we must instead calculate the three-body 
decay Z\ — > dqq. This decay is calculated in Ref. 26] via 7 and Z* exchange diagrams, but 
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Figure 3: Branching fraction of Z\ — > a + hadrons versus ■ 



neglecting interference terms. Here, we present the three-body width including interference: 



dT 



3 

12vr 3 



TTJ. _ v 1 



'I 



2g z Qeg v Re(G*G z )Re 



(1-Mfc) 2 ( 2 + Mfc) 



m - 



+ i 



(3.3) 



where the axino and quark masses have been neglected. In the above, Gz 



G 



v (1) 



gv = t£ — Q sin 2 6*iy and 5a 



16tt fa/N 

— T3 /2, where 



sin#iy, 



7 - IBS JJW cost, W, 9Z - S m6W cos0 w ' " T 

T3 is the weak isospin of the quark q. The above differential width is integrated over the 
range [i k : 4m 2 . /m 2 ~ — > 1. The quark mass acts as a regulator for the otherwise divergent 
photon-mediated contribution. 

The hadronic branching fraction Bh of Z\ is plotted in Fig. 



versus . The values 

^1 



of v\ , C a YY and f a /N cancel out in the branching fraction calculation, so the result is 



quite general. We see that for low values of < mz, Bh 



Once exceeds 
^1 



~ 0.02. 

mz, the branching fraction for Z\ — > dZ turns on and Bh increases asymptotically towards 
~ 0.175. Armed with the values of O^/i 2 , t(Z{), Bh and the constraints of Ref. [43], 
we are now ready to explore allowed regions of PQMSSM parameter space which might 
reconcile thermal leptogenesis with the gravitino problem and at the same time satisfy the 
BBN contraints for late decaying neutralinos. 



4. Allowed values of Tr in supersymmetric models 
4.1 Five benchmark models 

In this section, we discuss the sorts of models which might allow a reconciliation of thermal 
leptogenesis with the gravitino problem. Our starting point is to require models with a 
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rather heavy gravitino mg > 30 TeV, so as to avoid the gravitino BBN constraint. We 
also invoke mixed axion/axino dark matter with a light a (ma < 1 MeV) so as to avoid 
the gravitino-induced overclosure problem. 

Models with multi-TeV scale gravitinos are rather limited. In ordinary gravity- mediation 
(SUGRA) models, the gravitino mass arises from the superHiggs mechanism in a hidden 
sector of the model. The gravitino mass then sets the scale for the soft SUSY breaking 
terms. This latter condition applies to the scalar sector in simple SUGRA models, while 
gaugino masses require in addition stipulation of the gauge kinetic function. We will as- 
sume here that gauginos are quite light; if instead gaugino masses are in the multi-TeV 
range, then RGE effects drive the third generation scalar soft masses into the mufti- TeV 
range as well, thus engendering a conflict with naturalness. 

Two types of SUGRA models have multi-TeV (1st and 2nd generation) scalar masses, 
while the 3rd generation is around the TeV scale, as motivated by the hierarchy problem. 
The first is Yukawa- unified SUSY (YU), in which scalar masses are preferred to be in the 
mufti- TeV range at the GUT scale, while gauginos need to be quite light. In these models, 
the large, unified third generation Yukawa couplings drive the third generation soft terms 
down into the TeV range while first and second generation soft terms at Q = M wea k 
remain in the multi-TeV regime [45, 46, fi7[| . The SUSY particle mass spectrum for ji > 
is characterized as a radiatively driven inverted scalar mass hierarchy [48]. When requiring 
consistency with B-physics constraints, the lightest neutralino is a nearly pure bino state, 
while scalars are quite heavy thus suppressing the neutralino annihilation cross section. 
Generally, the neutralino relic density is large, of order Q,^ h 2, ~ 10 — 10 4 , far beyond 
the measured value. By invoking instead a light axino as LSP, the relic density of mixed 
axion/axino dark matter can be reconciled with observation [49, [50| . We present in Table [j] 



a Yukawa-unified benchmark model (model HSb from Ref. [52]) with m\§ = 10 TeV, which 
would be consitent with gravity mediation with a ~ 30 TeV gravitino mass. For this point, 
the lightest neutralino has mass ms = 49 GeV, and Vt^h 2 = 3195. 

The second type of gravity mediation model which would be consistent with ~ 30 TeV 
gravitinos is effective SUSY, or ESUSY |53[| . In Ref. [54], these models were explored with 
GUT scale soft SUSY breaking boundary conditions. Viable spectra with a weak scale 
inverted scalar mass hierarchy were found. For ESUSY models, first/second generation 
scalars could have mass mo (1,2) in the multi-TeV range at the GUT scale, while third 
generation scalar masses mo (3) are in the few TeV range. Upon RG evolution, first/second 
generation scalars remain in the multi-TeV range, while third generation scalar masses are 
suppressed by two-loop RGE terms, and are sub-TeV at the weak scale. Benchmark point 
BM2 is an example of an ESUSY scenario. 

If we proceed beyond gravity-mediation, then the class of models which necessarily 
supports multi-TeV gravitinos is anomaly- mediation (AMSB) |55|]. In these models, spar- 
ticle masses arise at the loop level via the superconformal anomaly. Sparticle masses are 
of order m ~ yg^m.3/2, so that in order to support TeV scale sparticle masses, a grav- 
itino mass of order 50 — 200 TeV is needed. In AMSB models, the lightest neutralino is 
nearly a pure wino state with a relic density usually well below the measured Odm^ 2 5 unless 
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Table 1: Masses and parameters in GeV units for the five benchmark points. BM 2-5 are computed 
with Isajet7.81 using m t = 173.1 GeV. BM1 uses Isajet 7.79 with m t = 172.6 to be consistent 
with previous work. 



nig > 1300 GeV. For benchmark model BM3, we select a gaugino AMSB (inoAMSB) point 
with 777,3/2 = 50 TeV. 4 It has been argued in Ref. [56] that in string models the scalar soft 
masses and trilinear terms are actually suppressed, while gaugino masses assume the usual 
AMSB form. These models, with too = Aq = at Mgut, avoid the problem of tachyonic 
scalars which occurs in traditional AMSB models; the scalar masses are uplifted via RG 
running during their trajectories from Mqut to M wea i-. While we do select the inoAMSB 
model as our benchmark point, very similar dark matter phenomenology occurs for mini- 



4 We use different notations for the gravitino mass scale m.3/2 and the physical gravitino mass mg, though 
« m 3/2 . 
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mal AMSB (mAMSB) or hypercharged AMSB [57], since the defining characteristics are a 
wino-like lightest neutralino [58]. 

The third class of models we examine also easily supports multi-TeV gravitinos. These 
are the mixed moduli- AMSB models ||59|], also known as mirage unification (MU). This class 
of models is inspired by the KKLT set-up of string models with flux compactifications and 
an uplifted scalar potential which can accommodate a positive cosmological constant p0|. 
While MU models require a multi-TeV gravitino mass, the lightest neutralino can easily 
remain bino-like, and can also have a very low relic abundance ^^Ji 2 at the 0.1 level or 
below. These models are stipulated by the parameters a, which governs how much gravity 
versus anomaly mediation occurs, along with 7713/2 and tan/3. One must also stipulate the 
matter and Higgs field modular weights n m and n# , which take on values of 0, \ or 1- The 
MU models are hard coded into the Isasugra spectrum generator. 

Benchmark model BM4 takes a = —1.6, 7713/2 = 30 TeV, tan {3 = 10 and n m = nu = 0. 
It has a bino-like lightest neutralino with = 146 GeV, but ^^^i 2 = 0.04 due to bino- 
wmo co- annihilation, or BWCA @. In BWCA, the gaugino masses M\ ~ —M2 at the 
weak scale. Since the gaugino masses have opposite signs, there is no mixing between bino 
and wino states, although they can be close in mass and can thus co-annihilate. If instead 
Mi ~ +M2 at the weak scale, then one obtains a Z\ of mixed bino-wino content (which 
also occurs in MU models). 

Benchmark model BM5 is also of the MU type, but with a = 6, 7773/2 = 30 TeV, tan f3 = 
10 and n m = ^, 77 h = 0. This model yields a bino-like Z\ with mass = 759 Gev, 
but Q^h 2 = 0.06 due to neutralino annihilation through the pseudoscalar ^4-resonance |62|| . 



In the following, we examine whether these scenarios are compatible with a Tr high 
enough to allow for thermal leptogenesis. To this aim we perform for each BM point a 
random scan over PQMSSM parameters in the range 

m- a e [10- 7 , 10" 1 ] GeV, 
f a /N € [10 8 , 10 14 ] GeV, (4.1) 
0i € [0, tt] . 

and calculate the value of Tr which is needed to enforce f2 a „ = 0.1123. As mentioned, a 
major constraint comes from the BBN bounds on the lifetime of the Z\. A digitized version 
of these BBN bounds is shown in Fig. ||, in the t{Z\) versus ^^Ji 2 plane. We also show 
the locus of benchmark points BM1-BM5 on the plot, along with the respective maximum 
allowed values of f a /N which are consistent with the BBN bounds. 

4.2 Gravity mediation: Yukawa unified SUSY 

The YU model point BM1 has = 49 GeV, tj^ = 1 and = 3195. In addition we 
take 777g = 30 TeV. Since ^j^/i 2 is so large, the Z\ lifetime is restricted to be < 0.03 sec by 
the BBN bounds from Ref. |43[| , see Fig. ||. The neutralino lifetime bound then translates 
into an upper bound on f a /N < 3 x 10 10 GeV. 
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Figure 4: BBN bounds on late-decaying neutral particles with = 0.1, digitized from Ref. [43 
in the t{Z\) versus ^^ 1 h 2 plane. We also show the locus of benchmark points BM1-BM5, along 
with their maximum allowed f a /N values. 




fJH (GeV) 

Figure 5: Scan over PQMSSM parameters for BM1 (YU model) plotted in the Tr vs. f a /N plane. 
The blue points respect the Z\ — > a + hadrons BBN bound, while the red points violate the BBN 
constraint. Points shown in light blue or light red have > 20% WDM or > 1% HDM, as discussed 
in the text. 



The results of the scan over PQMSSM parameters, Eq. Q4.1| ), are shown in Fig. [5] in 
the fa/N vs. Tr plane. For Tr values above the diagonal line labelled Tr = Td cp , the 
axinos would have been in thermal equilibrium in the early universe; all points lie below, 
so the expression for £l a ah 2 in Eq. (|2.7D is valid. In the figure, the [light and dark] blue dots 
denote points consistent with the neutralino BBN bound, while [light and dark] red dots 
denote points which violate BBN constraints. As shown in Fig. [l], the largest Tr values are 
mostly obtained for a light axino. Depending on its mass, the axino might constitute warm 
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(WDM) or hot (HDM) dark matter, which is severely constrained by the matter power 
spectrum and reionization |26|, [63|], see also [64, |S5]| . Since the bounds on the amount of 



HDM/WDM are model dependent |33|], we do not impose such constraints on our results. 
However, as a guidance, we indicate by lighter colors the points which have: 

• m- a < 100 keV and Q a /Q aa > 0.2 or 

• m a < 1 keV and Vt a /£l aa > 0.01, 

where Q a = ^ P + ^& + • Dark blue and dark red points thus have mostly CDM with 
at most 20% WDM and 1% HDM admixture. 5 

As can be seen in Fig. [|, for values of f a /N consistent with BBN bounds, only Tr 
values below 10 8 GeV are allowed. Such low values of Tr are insufficient to support thermal 
leptogenesis, but are sufficient to support non-thermal leptogenesis, which requires a more 



modest value of T R > 10 6 GeV g|, |§. Note also that these points (with T R ~ 10 6 - 
10 8 GeV) have a substantial fraction of WDM. In fact, as we will see later, for Z\ masses 
of about 50 GeV, as typical for YU scenarios with fi > 0, reconciling thermal leptogenesis 



with the gravitino problem requires a very low neutralino abundance, cf. Fig. 13. Such low 
abundances can indeed be achieved in a small region of the YU SUSY parameter space 
or-perhaps more easily-using either generational non-universality, or gaugino mass 



non- universality [51]. 



4.3 Gravity mediation: Effective SUSY (ESUSY) 

Another gravity-mediation model which requires multi-TeV scalars is effective SUSY. The 
ESUSY point BM2 has a mixed bino-higgsino NLSP with mass = 414 GeV. The Z^ 
and W\ are quite close in mass to the Z\, followed by i\ and 61 which are just 60% heavier. 
This leads to ^^Ji? ~ 0.04 due to simultaneous mixed bino-higgsino- wino enhanced anni- 
hilation, and also a contribution from stop and sbottom co-annihilation. The low ^^/i 2 
value allows Z\ lifetimes up to ~ 200 sec, corresponding to f a /N values as high as 10 13 
GeV (see Fig. |). 

Such high fa/N values suppress the thermal production of axino dark matter, while low 
values of 6i suppress the axion relic abundance. Scanning over PQMSSM parameters reveals 
that re-heat temperatures above 10 12 GeV can be generated while avoiding overproduction 
of dark matter and maintaining consistency with BBN bounds. Thus, ESUSY models with 
a low abundance of neutralinos, mixed axion/axino dark matter with a high PQ scale and 
low 6i, apparently can reconcile thermal leptogenesis with the gravitino problem, although 
most of the solutions for BM2 have a potentially dangerous fraction of HDM/WDM. 

We point out, however, that the BBN bounds will be less severe for similar scenarios 
with larger /j, (i.e. less Z\ higgsino admixture) but lighter t\ or b\, such that a low ^^/i 2 
arises from stop or sbottom co-annihilation; examples are discussed in |p5| . 

5 A rough estimate based on the neutrino mass limit |^5| from cosmological data, YJ rn„ < 0.41 to 0.44 eV, 
gives that up to 4-5% HDM contribution could be acceptable. 
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Figure 6: Scan over PQMSSM parameters for BM2 (ESUSY model) plotted in the T R vs. f a /N 
plane. Same color code as in Fig. |[ 

4.4 Anomaly mediation: gaugino AMSB 

Next, we consider a inoAMSB model, point BM3. In this case, the relic abundance is much 
lower, Q^h 2 = 0.016, thus much longer Z\ lifetimes of ~ 800 sec are allowed. Nevertheless, 

the Z\ lifetime is suppressed in this case by the tiny value of vjp ~ 0.01, so that for BM3, 
fa/N values only as high as 2 x 10 11 GeV are allowed. The results are shown in Fig. |^, 
where we again show BBN-allowed and BBN-forbidden model points in the f a /N vs. Tr 
plane. We see that just a few points barely exceed the rough requirement for thermal 
leptogenesis that Tr > 2 x 10 9 GeV. If we increase mg beyond 50 TeV, then the value 
of Q^h 2 increases, requiring shorter Z\ lifetimes, although the Z\ lifetime also decreases 

as ~ l/m~ . We also note that the light blue points with Tr > 2 x 10 9 GeV have axino 

1 

masses of a few times 10 -7 GeV. Overall, we conclude that AMSB models with a wino-like 
Z\ can just barely reconcile thermal leptogenesis with the gravitino problem, however this 
would require a considerable fraction of axino HDM. 

4.5 Mixed moduli/anomaly mediation with bino-wino co-annihilation 

Let us now move to mirage unification models, and the BM4 point with bino-wino co- 
annihilation. In this case, the neutralino mass is 146 GeV and its relic density f^/i 2 = 0.04, 
so that Z\ lifetimes of ~ 200 sec are allowed. Since the Z\ is nearly pure bino, its decay 
is unsuppressed by v^\ and we find in Fig. ^ that f a /N values as high as 10 13 GeV are 
allowed by BBN constraints. With such high f a /N values, the thermal production of 
axinos is suppressed, and Tr values over 10 12 GeV are allowed. Thus, these models are 
capable of reconciling thermal leptogenesis with the gravitino problem while avoiding BBN 
constraints. 

In Fig. we plot the scanned points for the BM4 point in the Q. L vs. Tr plane. (In 
fact the analogous plot for BM2 looks almost the same.) We see that the BBN-allowed 
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Figure 7: Scan over PQMSSM parameters for BM3 (gaugino AMSB model) plotted in the 
Tr vs. fa/N plane. Same color code as in Fig. ||. 
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Figure 8: Scan over PQMSSM parameters for BM4 (MM-AMSB with bino-wino co-annihilation) 
plotted in the Tr vs. f a /N plane. Same color code as in Fig. ||. 



blue points must have Q{ on the small side, certainly < 1 in order to avoid overproducing 
axions via vacuum mis-alignment, while maintaining Tr > 10 9 GeV. 



Figure 1C shows the scan results for BM4 in the rria vs. Tr plane. Here, we see that the 
points which are BBN-allowed, and also are consistent with reconciling thermal leptogenesis 
with the gravitino problem require rria < 100 keV. For values of mj lower than 100 keV, the 
axinos may start becoming warm rather than cold dark matter, and for values lower than 
1 keV they contribute to HDM. Thus, if CDM/WDM constraints are properly applied, we 
expect that a small region of parameter space will be consistent with thermal leptogenesis, 
as roughly indicated by the dark blue points in Figs. 
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Figure 9: Scan over PQMSSM parameters for BM4 plotted in the Tr vs. ftj plane. Same color 
code as in Fig. [sj. 




Figure 10: Scan over PQMSSM parameters for BM4 plotted in the Tr vs. mj plane. Same color 
code as in Fig. |5[ 



4.6 Mixed moduli/anomaly mediation with neutralino annihilation on the A- 
resonance 

Next, we turn to benchmark point BM5, a MU model with a bino-like neutralino with 
= 759 GeV and £lg h 2 = 0.06 due to neutralino annihilation through the ^4-resonance 
(itia = 1584 GeV, so that 2m^ ~ i^a)- The low value of Q^h 2 again allows for Z\ 
lifetimes as high as ~ 200 sec. But now, since is so large, f a /N values up to ~ 10 14 
GeV are allowed by Fig. ||. The scanned points are shown in Fig. [TI|. We see that all points 
pass the BBN constraints, allowing for Tr values in excess of 10 12 GeV. This model again 
easily reconciles thermal leptogenesis with the gravitino problem. 
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Figure 11: Scan over PQMSSM parameters for BM5 (MU with A-resonance annihilation) plotted 
in the Tr vs. f a /N plane. Same color code as in Fig. g. 
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Figure 12: Scan over PQMSSM parameters for BM5 MU model plotted in the Tr vs. mj plane. 
Same color code as in Fig. || 



As in the case of BM4, 9i values must be < 1, and for very high Tr > 10 11 GeV, 6{ < 
0.4. We conclude that small values of 0i are necessary to allow for thermal leptogenesis in 
SUSY models with mixed axion/axino dark matter. We also show in Fig. |l2]the scan result 
for BM5 in the mj vs. Tr plane. Here we see that consistency with thermal leptogenesis 
requires axino masses rria < 10 MeV. Thermally produced axinos with mass > 0.1 MeV 
should constitute cold dark matter, so these points would have a mix of cold axions plus 
cold thermally produced axinos. 
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5. Conclusions 



In this paper, we examined imparity conserving super symmetric models with a goal of 
reconciling thermal leptogenesis with the gravitino problem, via the postulation of mixed 
axion/axino dark matter. The mixed dark matter arises naturally from the Peccei-Quinn 
solution to the strong CP problem in super symmetric models. 

In order to reconciliate thermal leptogenesis with cosmological constraints, such as the 
dark matter relic abundance and the BBN bounds on late decaying particles, we conclude 
that the following conditions are necessary: 

• Tr > 2 x 10 9 GeV, to allow for efficient thermal production of right-handed neutrinos 
in the early universe; 

• m,Q > 30 TeV, to avoid BBN constraints on gravitino production in the early universe. 
Models with such a heavy gravitino are also favored in SUGRA models in that they 
set the scale for the scalar soft mass terms; if the scalar masses are sufficiently high, 
then they can suppress unwanted FCNC and CP violating processes and also proton 
decay via a decoupling solution |37j : 

• An axino LSP in keV to MeV mass range; this condition allows us to avoid overpro- 
duction of dark matter from axino thermal production in the early universe as well 
as from Z\ and G decays, since the matter density is then suppressed by the ratio 
nia/m^^ g. We also require a neutralino NLSP, which is common in supersymmetric 
models with heavy scalar masses; 

• fa/N > 10 12 GeV, to suppress thermal overproduction of axino dark matter. The 
high value of f a /N means the axion mass is likely to lie in the sub-micro-eV range; 

• @i ^ 1; to avoid overproduction of axions when f a /N > 10 12 GeV; 

• n^h 2 < 1, > 100 GeV and/or vjp ~ 1, to avoid BBN constraints on the late 
decaying neutralino. 

For the SUGRA case, we examined Yukawa-unified and effective SUSY models, since 
these can easily accommodate a 30 TeV gravitino mass. These models typically have too 
low an annihilation cross section for the neutralino NLSP in the early universe, which leads 
to conflicts with BBN constraints on late decaying neutral particles: in this case, hadronic 
neutralino decay via Z\ — > Z* /j — > aqq. 

The ESUSY model can more easily allow for low neutralino abundances via stop, 
sbottom, stau or higgsino co-annihilation. In these models, requiring the sum of four 
production mechanisms for mixed axion/axino DM to equal the measured abundance can 
allow for Tr > 10 10 GeV, thus reconciling thermal leptogenesis with the gravitino problem. 

In AMSB models, one naturally has a gravitino mass in the 50-100 TeV range and small 
Qg^h 2 , but with a wino-like neutralino. The decay rate of the wino-like Z\ is suppressed 

(1) 2 ~ 

by the mixing factor v\ , leading to long-lived Zis, and likely conflicts with BBN. 
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We also examined models with mixed moduli-AMSB (mirage unification) soft terms. 
These models allow for 30-100 TeV gravitinos, but with a bino-like neutralino, so its lifetime 
is typically < 100 sec. We examined two cases: bino-wino co- annihilation and ^4-resonance 
annihilation. Both cases easily allow Tr to reach over 10 12 GeV, thus easily reconciling 
thermal leptogenesis with the gravitino problem, while respecting BBN constraints on 
long-lived neutralinos. 

These findings are summarized in a model-independent way in Fig. 13, which shows 
PQMSSM scan points with T R > 2 x 10 9 GeV in the % x /i 2 vs. f a /N plane for = 

50 and 500 GeV. 6 In both cases, = 1; recall that the neutralino lifetime scales as 
[v^ /(fa/N)] 2 . It is interesting to note that the reconciliation of thermal leptogenesis 
and the gravitino problem, in the framework used for our analysis, spans a wide range of 
Qg^h 2 , from Sl^/i 2 ~ 1 down to very low values. In particular for a light neutralino NLSP, 
very low Q^h 2 is required. Models with neutralino DM, on the other hand, would have 
Qj^/i 2 ~ 0.1, while models with a light axino DM and non-thermal leptogenesis would 
prefer ^^_h 2 > 100. Therefore distinct regions of the (PQ)MSSM parameter space are 
prefered depending on which DM and baryogenesis solutions are chosen. 

Experimental consequences of this scenario to reconcile thermal leptogenesis with the 
gravitino problem in super symmetric models include 1. a discovery at LHC of any of the 
models discussed here (or others) which support a gravitino mass in excess of 30 TeV, 2. the 
inferred apparent relic abundance of neutralinos is typically Qgji 2 < 1, 3. null results from 
direct or indirect WIMP searches, and 4. a positive signal for the QCD axion at sub-/ieV 
levels at ADMX [66] or other axion detection experiments. 
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